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Abstract. Wc report the results of simultaneous obser- 
vations of the Vela pulsar in X-rays and radio from the 
RXTE satellite and the Mount Pleasant Radio Observa- 
tory in Tasmania. We sought correlations between the 
Vela's X-ray and radio flux densities and radio arrival 
times on a pulse by pulse basis. We found significantly 
higher flux density in Vela's main X-ray peak during ra- 
dio pulses that arrived early. This excess flux shifts to 
the 'trough' following the 2nd X-ray peak during radio 
pulses that arrive later. We suggest that the mechanism 
producing the radio pulses is intimately connected to the 
mechanism producing X-rays. Current models using reso- 
nant absorption in the outer magnctosphcre as a cause of 
the radio emission, and less directly of the X-ray emission, 
are explored as a possible explanation for the correlation. 



1. Introduction 

The complexity of Vela's spectrum allows for the possi- 



bility of both polar cap (Daugherty & Harding 1996) and 



outer-gap (Cheng, Ruderman, & Zhang 2000) models of 



emission, but only a subset of the models suggest a connec- 
tion between the radio and X-ray emission. For example 



Harding et al.(2002) suggest that scattering of soft X-ray 



photons in a polar cap model could explain the observed 
optical spectrum of Vela and also provide for the align- 
ment of the radio peak with a very soft X-ray peak, which 
they observe. The connection between X-ray and radio 
emission that wc offer in this article favors models such as 
this, which connect far reaching ends of the electromag- 
netic spectrum. 

Additional experiments linking pulsar emission in 
the radio and X-ray regimes have been performed 
by Cusumano ct al. (2003) and Vivekanand (2001). 



Cusumano et al. show that in PSR B1937+21 there is 
close phase alignment between X-ray pulses and giant ra- 
dio pulses, suggesting a correlation in their emission re- 
gions. The work by Vivekanand, on the other hand, ex- 
plains the connection between the Crab pulsar's radio and 
X-ray flux variations as a result of temporal variations in 
the number of "basic emitters" or in the basic emitters' 
degree of coherence. 

Giant radio pulses have not been shown to exist in 
Vela, but Johnston et al. (2001, hereafter J01) discovered 
'giant micropulses' in the Vela pulsar, which have flux den- 
sities no more than ten times the mean flux density and 
have a typical pulse width of ~ 400/is. These micropulses 
may explain the results of Krishnamohan & Downs (1983, 
hereafter KD83) who found that the strongest radio pulses 
arrive earlier than the averaged profile. 

By doing a pulsc-by-pulse analysis of the Vela pulsar 
in X-ray and radio wavelengths, we will show in this paper 
that the emission mechanisms creating the Vela pulsar's 
X-ray and radio flux densities must be related. We will 
discuss the X-ray and radio observations in |j2j our anal- 
ysis in |j3l the effects of scintillation in Sj4j a discussion of 
interpretations in [}5l theoretical predictions as they relate 
to our results in Sj6l and finally our conclusions and related 
future work in fJ71 

2. Observations 

Our data consist of 74 hours of simultaneous radio and 
X-ray observations taken over three months at the Mount 
Pleasant Radio Observatory in Tasmania and with the 
RXTE satellite. 

The radio data were acquired during 12 separate obser- 
vations using the 26m antenna at a frequency of 990.025 
MHz between 30 April and 23 August, 1998. All individ- 
ual pulses from Vela are detectable, and the pulse height, 
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integrated area, and central time of arrival (for the solar- 
system barycenter) were calculated from cross-correlation 
with a high signal to noise template in the usual fashion. 

The X-ray data were taken during the same three 
months, yielding 265 ks of usable simultaneous obser- 
vation. For the purposes of this project, only top-layer 
data from RXTE's Proportional Counter Units (PCUs) in 
Good Xenon mode in the energy range of 2-16 keV were 
used. Other filtering parameters included were standard 
RXTE criteria: elevation was greater than 10 degrees, off- 
set was less than 0.02 degrees, the data were taken with at 
least 3 PCUs on, time since SAA was at least 30 minutes, 
and electronO was less than 0.105. 



3. Analysis 

We filtered the X-ray photon arrival times and trans- 
formed them to the Solar System Barycenter (SSB) us- 



ing the standard FTOOLS (Blackburn 1995) package. We 



calculated the pulsar phase at the time of each X-ray pho- 
ton, and matched each X-ray photon with the radio pulse 
that arrived at the SSB at the same time. We then com- 
pared pulse profiles for X-rays segregated according to the 
arrival time of the radio pulse. 

We filtered the X-ray photon arrival times and trans- 
formed them to the Solar System Barycenter (SSB) using 
the standard FTOOLS |(Blackburn 1995)| package. We cal- 
culated the pulsar phase at the time of each X-ray photon, 
using the radio pulsar-timing program TEMPO 1 , and the 
ephemeris downloadable from Princeton UniversitjQ. We 
matched each X-ray photon with the radio pulse that ar- 
rived at the SSB during the same turn of the pulsar. The 
precise period of time associated with each radio pulse was 
given by our best model for the pulse arrival time, ±0.5 x 
the instantaneous pulse period calculated via the model. 
Photons arriving on the borderline were associated with 
the earlier pulse. We then compared pulse profiles for X- 
rays segregated according to the corresponding radio pulse 
arrival time. 

Radio pulses arrive at a range of times around the 
predicted arrival time, as KD83 found. The histogram of 
residual arrival times for radio pulses, relative to the pre- 
diction of our best model, is shown in FigureQ] We divided 
all of the pulses into four quartilcs, by the residual phase 
of the radio pulse, with equal numbers of pulses in each 
quartile. Figure Q] shows our division of the residual phase 
of the radio pulse into the quartiles. We called the quartiles 
early, medium-early, medium-late, and late, according to 
the arrival time of the radio pulse. 

We formed an X-ray pulse profile for each of the quar- 
tiles of radio-pulse arrival times, from the X-ray photons 
associated with each. Figure [2] shows the four resulting 
X-ray profiles with errors determined by counting statis- 
tics. The X-ray profiles are significantly different; the pulse 



changes in shape among the quartiles. In particular it ap- 
pears that the first X-ray pulse is sharper and stronger 
during "early" radio pulses. Also, it appears that the 
trough after the 2nd X-ray peak around phase 0.4 (here- 
after called simply 'the trough') gets "filled in" in late 
quartiles. 

In order to graphically analyze the changes in pulse 
shape, we plot the profiles in different quartiles against 
one another in Figure [3l The solid diagonal line shows 
where the points would lie if the profiles were identical. 
As the figure shows, Bin 2 (the tallest bin in the 1st X-ray 
peak) is significantly stronger in the early quartile than 
in the late quartile. In addition, although less noticeable, 
Bins 12, 13 and 1 (the trough) are significantly weaker in 
the early quartile than in the late. 

Overall, the data suggest that sharper X-ray pulses are 
associated with earlier radio pulses. In Figure |4] we plot 
the height of the peak in Bin 2, and of the trough in Bins 
12 and 13, against the lateness of the quartile. The figure 
shows a systematic trend of a weaker peak toward late 
arrivals. The increase in photon rate in the trough makes 
up for much of the loss of X-ray photons in the peak. 
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See http:/ /pulsar. princeton.edu/tempo 



Fig. 1. The total number of radio pulses vs phase relative 
to a predictive long-term timing model. The dotted lines 
show the boundaries of the four bins that were used to 
make the four profiles shown in Figure [2] 



4. Scintillation 

In contrast to effects intrinsic to the pulsar, scintillation 
is unlikely to produce the observed association, because 
it does not affect X-rays; scintillation might erase such 
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Fig. 2. Full-period X-ray profiles for photons detected 
during radio pulse arrival times falling in the four quartile 
bins shown in Figure [T] The radio peak falls at phase 0.5. 
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Fig. 3. Measured counts in the late radio phase quartile 
vs counts in the early radio phase quartile. X-ray phase 
bins are numbered sequentially beginning with 1. 



a correlation but it cannot introduce it. Nonetheless, we 
used a number of techniques to ensure that our result 
was not produced by scintillation including subtracting 



Fig. 4. The height of Bin 2, and of the dip in Bins 12 and 
13, vs the lateness of the quartile. 



the average residual from each 5-minute span of data, 
and defining the four quartiles based on both 5-minutc 
and 1-hour spans of data. More details can be found in 
Lommen et al.(2006) In summary, none of these tech- 



niques changed our results significantly. 



5. Discussion 

J01 showed that the giant micropulse emission occurs 
about 1 ms before the radio peak, so it is realistic to con- 
sider the possibility that the giant micropulse emission 
is primarily responsible for the early arrival of the radio 
pulse. The hypothesis is consistent with the timescale of 
the "early-ness" : the full- width at half- maximum of the 
histogram of arrival time residuals shown in Figure [1] is 
about 0.5 ms. However, the more detailed analysis per- 



formed in Lommen et al.(2006) shows that the effect can 



not be explained by the occurrence of giant micropulses. 
6. Theoretical Predictions 

The results described above imply a connection between 
the radio and X-ray emission mechanisms for Vela that 
is not consistent with outer gap models. In these models, 
the high energy emission results from a gap connection to 
the pole opposite from that producing the radio emission. 
It is not clear how a correlation could exist between the 
radio and high energy regimes in these models. Pctrova 
(2003, and references therein), however, offers a detailed 
model that could explain the correlation. She states that 
the non-thermal optical emission of rotation-powered pul- 
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sars should be considered as part of the broadband high- 
energy emission. Her models suggest that resonant absorp- 
tion of radio emission from the outer magnetosphere leads 
to an increase in the pitch angles and momenta of the 
secondary pairs, which then leads to optical and higher 
energy emission by spontaneous synchrotron radiation. 
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7. Conclusions and Future Work 

We have significant evidence linking features of Vela's X- 
ray emission with features of its radio emission. We find 
that X-ray pulses associated with early radio pulses show 
stronger emission at the main X-ray peak which is the 
sharper of the two. Similarly X-ray pulses associated with 
later radio pulses show stronger emission at the trough 
following the 2nd X-ray peak. We conclude that there is 
a close relationship between X-ray and radio emission in 
the Vela pulsar. 

We plan to further characterize the relationship be- 
tween the radio and high energy emissions of pulsars to 
identify their origins and constrain magnetospheric mod- 
els. In particular, we will explore the dependence of radio- 
to- X-ray correlations on the radio frequency and polariza- 
tion properties of individual Vela pulses, both of which 
carry information about emission altitudes. Similar ob- 
servations of other pulsars also promise useful insights 
as probes of different magnetic field strengths and emis- 
sion/viewing geometries. 
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